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Abstract 

The correlation between cosmic shear as measured by the image distortion of high-redshift 
galaxies and the number counts of foreground galaxies is calculated. For a given power 
spectrum of the cosmic density fluctuations, this correlation is proportional to the bias 
factor, which can thus directly be measured. In addition, this correlation provides a first- 
order measure of cosmic shear and is therefore easier to observe than quadratic measures 
hitherto proposed. Analytic approximations show that the expected signal-to-noise ratio 
of the correlation is large, so that a significant detection is possible with a moderate 
amount of data; in particular, it is predicted that the ongoing ESO Imaging Survey 
(EIS) will be able to detect this correlation on scales of ~ 10' at a 3-cr level, and at with 
higher significance on smaller angular scales. 

1 Introduction 

The distortion of high-redshift galaxy images by the (tidal) gravitational field of inter- 
vening matter inhomogeneities (often called 'cosmic shear') can be used to study the 
intervening mass distribution. In particular, if the large-scale structure of the (dark) 
matter is considered, the observable image distortions constrain the statistical properties 
of the cosmic matter distribution. This method of determining the power spectrum of cos- 
mic density fluctuations has been investigated recently in considerable detail (Blandford 
et al. 1991; Miralda-Escude 1991; Kaiser 1992, 1996, hereafter K92, K96, respectively; 
Villumsen 1996; van Waerbeke, Bernardeau & Mellier 1997; Jain & Seljak 1997; Sanz, 
Martinez-Gonzalez & Benitez 1997; Schneider et al. 1997a, hereafter SvWJK). A first 
significant detection of cosmic shear has been reported in Schneider et al. (1997b). 

Cosmic shear probes the statistical properties of the projected density fluctuations, 
where the projection takes into account the redshift distribution of the source galaxies 
and geometric factors. Similarily, the surface number density of galaxies is obtained by 
a projection of the three-dimensional galaxy distribution. Provided galaxies trace the 
underlying dark matter distribution, these two projected fields are correlated. 

In this letter, this correlation is investigated, aiming at a method to constrain the 
bias factor b which relates the number density fluctuations of galaxies to those of the 



underlying dark matter distribution. 1 A brief summary of the aperture mass statistics 
as a measure for cosmic shear is presented in Sect. 2, and an analogous statistics is 
introduced for the galaxy number counts. The general expression for the correlation 
between these two measures is derived in Sect. 3, and practical estimators are considered 
in Sect. 4. Focusing on large angular scales, linear theory presents a useful approximation 
for the growth of cosmic density fluctuations; in Sect. 5, the general expressions will be 
evaluated in this approximation. In particular, it is shown that the signal-to-noise ratio 
for the correlation coefficient is of order unity even in a single field. Therefore, this 
correlation should be easily detectable in currently conducted wide-field surveys, such 
as the EIS (Renzini & da Costa 1997). In particular, the significant verification of this 
correlation is probably the easiest way to detect cosmic shear. 

2 Aperture mass and number counts 

In this section we briefly summarize the M ap -statistics for cosmic shear, and introduce a 
similar statistics for the number counts of (foreground) galaxies, following the notation 
of SvWJK. 

Light propagation through a slightly inhomogeneous Universe can be described by 
an equivalent single-plane gravitational lens equation, to first order in the Newtonian 
gravitational potential (see SvWJK for a detailed discussion of this point) . For sources 
with a redshift probability density p z (z) dz = p w (w) dw, where w is the comoving distance 
out to redshift z, the dimensionless surface mass density at angular position 6 of this 
single-plane lens is 

= o — ) Qa / dw 9{w)fK{w) — , (1) 

2 V c J J a{w) 

where 5 is the density contrast, a = (1 + z) -1 the cosmic expansion factor, f2 d the density 
parameter in dust at present, g(w) := f w H dw' p w (w')fK(w' — w) / fx{w') is the source- 
averaged distance ratio D ds /D s for a density fluctuation at distance w, fx{w) is the 
comoving angular diameter distiance to comoving distance w, and wh is the comoving 
distance to the horizon. 

The Jacobi matrix, which describes the locally linearized lens mapping, reads 
= $ij ~ VVj'W) where indices preceded by a comma denote partial derivatives 
with respect to the components of 6, and ifj is related to k via the Poisson equation, 
V 2- i/> = 2k. The two components of the shear, here written in complex notation, are 
derived from the traceless part of A, 'j(O) = (V'.n ~ V',22) /2 + iV 7 , 12 - The shear therefore 
describes the tidal part of the deflection potential which causes the distortion of images. 

1 A different method to obtain the bias factor directly lensing uses the magnification effect which 
causes a correlation between foreground galaxies and high- redshift QSOs (e.g., Bartelmann 1995, 
Sanz et al. 1997, Dolag & Bartelmann 1997) and changes the angular two-point correlation function 
of high-redshift galaxies (e.g., Moessner, Jain & Villumsen 1997). 
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Provided the density contrast S is a homogeneous and isotropic random field, so is 
the projected density k. The power spectrum P K (s) of k is related to the power spectrum 

— * 

P(k) of the density fluctuations 5 through 



4 V c y J a 2 (w) \fK{w) J 

see K92 and K96 for a derivation of (2). The second argument of P indicates that the 
power spectrum evolves with redshift. Several sample power spectra P K are plotted in 
Fig. 1 of SvWJK. 

In SvWJK, the aperture mass 

M ap (6):= J d 2 tf£/(|tf|)«(tf) (3) 

was introduced as a statistics for measuring cosmic shear. Similar quantities had previ- 
ously been considered in somewhat different contexts (e.g., Fahlman et al. 1994; Kaiser 
1995; Kaiser et al. 1994; Schneider 1996). Here, 11(1}) is a compensated filter function, 
i.e., J Q d# = 0, which vanishes for $ > 9. The definition (3) is particularly useful 

since M ap can directly be expressed in terms of the shear, 

M ap (#) = J d 2 tfQ(|tf|) 7 t(tf), (4) 

where Q(#) = (2/# 2 ) Jq d$' §' U ($') — U ($), and the tangential component of the shear 
at a position -d = (#cos</?, #sin</?) is 7t(i?) = —Tie (7(1?) e -2 "''). Hence, on the one hand, 
-^ap yields a spatially filtered version of the projected density field, and on the other 
hand, it can be expressed simply in terms of the shear. Since in the weak lensing regime, 
the observed galaxy ellipticities provide an unbiased estimate of the local shear, M ap is 
directly related to observables. The dispersion of M ap is related to the power spectrum 
PM by 



poo 

M 2 ap (6)) = 2n / dssP K (s)[l(s9)} 
Jo 



(5) 



where I(rj) = dx x u(x) 3 (r]x)] here, we have written U(fi) = u(i}/8)/6 2 , and J n (x) 
denotes the Bessel function of first kind. Hence, M ap (6>) provides a filtered version of the 
projected power spectrum, and the width of the filter, here expressed by I 2 , depends on 
the choice of u. 

Provided that galaxies are biased tracers of the underlying (dark) matter distribu- 
tion, the expected number density of galaxies in the direction 6 is given by (cf. Bartel- 
mann 1995; Dolag & Bartelmann 1997; Sanz et al. 1997) 



N{0) = N 



1 + 6 J dw p g (w) 6 (f K (w)0,w) 



(6) 



where p g (w) is the probability distribution of the galaxies in comoving distance (or, 
equivalently, redshift), which depends on the selection criteria of the galaxy sample (such 
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as limiting magnitude, color, etc.), N is the mean number density, and b is the average 
bias factor for this galaxy sample. In analogy to the aperture mass, we define the aperture 
number counts 

N{9) = J d 2 $U(\D\)N(tf) , (7) 
with the same function U as in (3). 



3 Density — shear correlations 

The correlation between M ap (6>) and Af(9) is measured by 

C(9):=(M &p (9)M(9)) , (8) 

where the angular brackets denote the ensemble average. Inserting the explicit expres- 
sions for M ap and A/", and using the same method as in K96 to calculate the resulting 
projection of the correlator (55), a few manipulations similar to those in SvWJK yield 

m=i J^)\ db N l im 2MjM f^.pf' \ m] >. (3 . 4) 
V c / J a(w)f K (w)J \Jk{w) J 

Hence, the correlation C depends on the cosmological model, the redshift distributions of 
the galaxies which are used to estimate the shear (which we shall call 'background' galax- 
ies in the following, though this should not imply that all these galaxies are lying behind 
those from which J\f is measured) and those with which M is estimated ('foreground' 
galaxies). In particular, C is proportional to the bias factor b. In Sect. 5 below, we shall 
calculate C for a simple model of the power spectrum of cosmic density fluctuations, but 
first we turn to practical estimators of C. 



4 Practical estimators 

Assume that in a circular aperture of radius 9 there are N h galaxies at positions #j 
whose ellipticity is measured and which are thus used to estimate the shear, and that Nf 
galaxies with positions cfij are used for measuring M. An estimator for M ap is 

irf) 2 Nh 

^ap(60 = ^£Q(l*l)eti, (10) 
b i=l 

where e t i is the tangential component of the image ellipticity, defined in analogy to 7 t . In 
the limit of weak lensing, to be considered here, the relation e = e^ s ^ + 7 between image 
ellipticity e and intrinsic source ellipticity e^ s ^ holds, so that is an unbiased estimate of 
7(i?i) because the intrinsic orientations of the source galaxies are assumed to be random. 
An estimator for Af is 
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Nf 

Af(6) = J £u(\c Pj \), (11) 

k=i 

so that an estimator for C reads C(9) = M ap (6>) N(6). To obtain the expecta- 
tion value E((7), several averages have to be taken. First, one has to average over 
the intrinsic source ellipticity distribution. Denoting this operator by A, one finds 
k(c(9) S j = A ^M ap (6>)j Af(9), since Sf is unaffected by A. The next average has 
to be taken over the galaxy positions. The corresponding operator P factorizes into 
two operators Pi and P 2 , which read Pi = / I^S p 2 = UkLi I d VfcPOfc) 

where it has been ussumed that the 'background' galaxies are distributed randomly 
in angle, and the probability density p(iPk) f° r finding a 'foreground' galaxy at <p k is 
P(Vfc) = N((p k )[f d 2 f N^cp)] -1 ~ N((p k )/N{. Thus, we neglect deviations of the to- 
tal number of 'foreground' galaxies in the circular aperture from the expected num- 
ber; these deviations are of minor importance only, provided Nf 1. Then, by per- 
forming both averages, one finds that the expectation value of C(9) is indeed C(9), 
E(C(9)) ee <P(A(C(0)))) = C{9). 

We next consider the dispersion of C for the case that the 'foreground' galaxies 
are unrelated to the matter distribution which distorts the background galaxies, or in 
other words, that M ap and Af are uncorrelated. In that case, the expectation value of C 
vanishes, and 

ol := E (C 2 ) = E (M a 2 p ) E (>) . (12) 
For the first of these factors, one finds with the same methods as used in SvWJK that 

E(M a 2 p ) = (Mi p (9)) + ^ , (13) 

where G = n9 2 J d 2 i} (5 2 (|i9|), cr e is the dispersion of the intrinsic galaxy ellipticities, and 
a subdominant 'shot-noise' term has been dropped. For the second term in (4.9), one 
finds 

E(Af 2 ) = (Af 2 (9)) + 1g , (14) 

where G = n9 2 J d 2 $ U 2 (\il)\). The dispersion of M can be obtained either directly from 
observations, or can be calculated following the biasing hypothesis. With steps very 
similar to those used for deriving (9), one obtains 

(A/- 2 (#)) = 27r& 2 iV 2 / dw^\ [ dssp(-^—-,w]l 2 (9s). (15) 
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5 Analytical estimates 



This section presents analytical estimates for C and the corresponding signal-to-noise 
ratio, which are obtained after several simplifications. First, we shall assume an 
Einstein-de Sitter (EdS) cosmological model. Second, only angular scales larger than 
~ lOarcminutes are considered; for an estimate on such large angular scales, the 
power spectrum of the density fluctuations can be assumed to evolve linearly, so that 
P(k, w) = a 2 (w) P(k, 0) = a 2 (w) Po(k). Third, since by a convenient choice of the func- 
tion £/"(#), the resulting filter function I 2 is quite narrow (see Fig. 2 of SVwJK), only a 
small range of wavenumbers contribute to the integrals over the power spectrum; hence, 
we can locally approximate the power spectrum by a power law in k. And finally, we 
shall choose the redshift distributions of sources and 'foreground' galaxies to be very 
localized. 

Thus, let Po{k) = Ak n , with a slope n, and an amplitude which shall be de- 
termined from the rms fluctuations as in a sphere of radius R = 8/i _1 Mpc. Then, 
A = 2n 2 (i 1 aiR ( - 3+n \ with Ci = 9 J °° dx x n " 4 [sin(x) - xcos(x)] 2 . For EdS, one finds 
Jk{w) = w = (2c/Ho)[l — (1 + z) -1 / 2 ], wh = 2c/Hq. We assume that the background 
sources are all at the same redshift z 8 , so that p z (z) = 5n(z — z s ), and g(w) = (w s — w)/w s 
for w < w s and zero otherwise, and w s is the co moving distance out to redshift z s . Sim- 
ilarily, we assume that the foreground galaxies are well localized around the redshift z s , 
corresponding to comoving distance w g . Except for the calculation of (A/" 2 ) below, we 
shall approximate p g by a delta 'function'. 

The compensated filter function U (i?) will be the same as in SvWJK, namely 2 u(x) = 
(9/7r)(l -x 2 )(1/3-x 2 ) for which the corresponding function Q(-&) = q($/9)/9 2 is q{x) = 
(6/tt) x 2 (1 - x 2 ). Then, I(rj) = (12/tt) (J 4 (??)/?7 2 ), and G = G = 6/5. 

With these assumptions and simplifications, we now estimate C{9) according to (9): 
Defining dimensionless comoving distances by w = wH /c, (9) becomes 

C{9) = ^^ bNa 2 aK) g) 9-^ , (16) 

Ci w s Wg n ' 

where r = RH /c = (8/3) x 10~ 3 and (2 = J °° dx x 1+n I 2 (x). We shall write (16) in the 
following convenient form, 



C0» = 1.822 x 10- CM b ft *l ttK> (^)" < " +2 ' , ("I 



where a good approximation of £(n) is ((n) wl + 1.04(n + 3/2) + 0.275(n + 3/2) 2 , which 
is accurate to better than 0.2% for n G [—2, —1]. 

Next we calculate (A/" 2 ) from (15) with the assumptions listed above. However, the 
occurrence of the factor pl{w) in the integrand of (15) shows that we cannot simply use 



This choice corresponds to I = 1 in SvWJK 
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a delta-function distribution of the redshifts of the 'foreground' galaxies. The reason 
for this problem lies is the derivation of (15): it is valid only if the functions which 
are projected vary on a much larger scale than the largest wavelength on which the 
power spectrum has appreciable amplitude. This condition is obviously violated if the 
'foreground' galaxies are very sharply localized. Therefore, we assume that p g is constant 
in an interval of width uw g around w g , with v <C 1. Then, 



2/ \ / a \ -("+2) ( 18 ) 

= 1.215 x 10- «„)^^I^) (^) 

We compare the two terms in (14) to see which one dominates. Therefore, consider the 
ratio 

1 {N\d)) 7T0 2 _ i _ _ 6 2 al a> g ) / N \ ( 6 



Ai ' NG 



We thus see that Ai is typically smaller than unity in a situation when v <C 1 and w g < 1 
(note that i& g = 1/3 corresponds to z g = 0.44). 

For the same model, (M 2 p (0)) can be calculated. One finds: 

v -(2+n) 

<M a 2 p (#)> = 2.733 x lO" 3 C(n) a 2 Z(n) ^-j , (20) 

where Z(n) = J* dx (1 — x) 2 x~ n = 2/[(l — n)(2 — n)(3 — n)]. Comparing the two terms 
in (13), we define the ratio 

1_ _ (MjjM) 2N b 

A * °> G (21) 

SV ; 8 s V 0.05 yV 107 V20arcmin- 2 ; V0.27 

where Z(— 1.5) was taken as a fiducial value. We therefore conclude that the first term 
in (13) dominates in all cases of interest here (i.e., for angular scales larger than ~ 10'), 
A2 <C 1. With (18) and (20), the signal-to-noise ratio becomes 

l-n 

C / «v\ f w<?\ 2 / v 1 



°o V 1 w s )\w s ) V Z(n) v /(i + A 1 )(l + A 2 ) ' (22) 

with Ai < 1, A 2 1 in typical situations. This ratio is indeed encouragingly large: 
Consider, for example, foreground and background galaxies with distance ratio w g /w s = 
1/2; then, for a width parameter v = 0.2 and spectral index n = — 1.5, the signal-to- noise 
ratio is 0.42 times the A-dependent terms. Note that the dependence of the signal-to- 
noise ratio on the number density of foreground and background galaxies enters only 
through the A;-factors; as long as A; <C 1, C/a is independent of these densitites. 
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6 Discussion 

In this letter, a statistical measure for the correlation C of cosmic shear with the number 
density of 'foreground' galaxies was defined and calculated in terms of the power spec- 
trum of cosmic density fluctuations. A practical unbiased estimator for this correlation 
was defined, and its dispersion calculated. On large angular scales, linear theory yields 
an accurate estimate for the power spectrum of cosmic density fluctuations; in the frame- 
work of this approximation, C and the corresponding signal-to-noise ratio were evaluated 
explicitly. 

A measurement of C would yield, for each cosmological model and initial power 
spectrum P{k), a direct estimate of the bias factor b averaged over angular scale 9. The 
method allows to probe the scale and redshift dependence of the bias parameter. On a 
short term, the detection of cosmic shear via the correlation C is perhaps more useful: 
assuming the validity of the biasing hypothesis, C is a first-order measure of the cosmic 
shear. The large signal-to-noise ratio (22) per single field shows that it should be much 
easier to get a significant detection of C than for the previously proposed quadratic 
estimators of the shear. 

Taking the EIS as an example, with 'foreground' galaxies chosen to have I < 21, and 
'background' galaxies with 22 < I < 23.5, the number densities will be approximately 
3 and 7 per arcmin 2 , and the characteristic redshifts of the two galaxy populations will 
be ~ 0.3 and ~ 0.8, respectively (e.g., Lilly et al. 1995). Hence, even in this case the Aj 
are smaller than unity, and the signal-to-noise ratio per field will be of order 0.3. Thus, 
with ~ 100 fields taken from the EIS, a 3-a detection of the correlation on angular scales 
larger than ~ 10' should be possible, provided the data are of sufficient image quality. 

In a future publication, C will be calculated on smaller angular scales, using the 
fully non-linear evolution of the power spectrum (e.g., Peacock & Dodds 1996), and 
for different cosmological models. It is expected that the signal-to-noise ratio is not 
strongly affected by the non-linear evolution, and that an accurate measurement of C(6) 
on small angular scales (~ few arcminutes) will be possible with a moderate amount of 
high-quality image data. 

I would like to thank L. da Costa, E. Martinez-Gonzalez & L. van Waerbeke for 
fruitful discussions, and M. Bartelmann & S. Seitz in addition for helpful comments on 
the manuscript. This work was supported by the "Sonderforschungsbereich 375-95 fiir 
Astro-Teilchenphysik" der Deutschen Forschungsgemeinschaft. 
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